Lineage-specific Gene Duplications Research Articles

Copy number and selection of MHC genes in ruminants are related to habitat, average life span and diet

The ruminants, as the main group of livestock, have been extensively studied in terms of their physiology, endocrinology, biochemistry, genetics, and nutrition. Despite the wide geographic distribution and habitat diversity of animals in this group, their ecology and evolution remain poorly understood. In this study, we analyzed the gene copy number, selection, and ecological and evolutionary processes that have affected the evolution of major histocompatibility complex (MHC) genes across ruminant lineages based on available genomic data. The 51 species analyzed represented all six families of ruminants. Our finding indicated that the architecture of the MHC region is conserved in ruminants, but with variable copy numbers of MHC-I, MHC-IIA, and MHC-IIB genes. No lineage-specific gene duplication was observed in the MHC genes. The phylogenetic generalized least squares regression (PGLS) model revealed association between ecological and biological factors (habitat and lifespan) and gene duplication in DQA and DQB, but not in DRB. The selection pressure of DQA and DQB were related with lifespan, diet, and the ratio of genetic repeat elements. These results suggest that the MHC evolution in ruminants, including copy number and selection, has been influenced by genetic repeat elements, pathogen exposure risk, and intrinsic cost of possessing multiple MHC genes.

Identification and characterization of the cupin_1 domain-containing proteins in ma bamboo (Dendrocalamus latiflorus) and their potential role in rhizome sprouting.

Cupin_1 domain-containing protein (CDP) family, which is a member of the cupin superfamily with the most diverse functions in plants, has been found to be involved in hormone pathways that are closely related to rhizome sprouting (RS), a vital form of asexual reproduction in plants. Ma bamboo is a typical clumping bamboo, which mainly reproduces by RS. In this study, we identified and characterized 53 Dendrocalamus latiflorus CDP genes and divided them into seven subfamilies. Comparing the genetic structures among subfamilies showed a relatively conserved gene structure within each subfamily, and the number of cupin_1 domains affected the conservation among D. latiflorus CDP genes. Gene collinearity results showed that segmental duplication and tandem duplication both contributed to the expansion of D. latiflorus CDP genes, and lineage-specific gene duplication was an important factor influencing the evolution of CDP genes. Expression patterns showed that CDP genes generally had higher expression levels in germinating underground buds, indicating that they might play important roles in promoting shoot sprouting. Transcription factor binding site prediction and co-expression network analysis indicated that D. latiflorus CDPs were regulated by a large number of transcription factors, and collectively participated in rhizome buds and shoot development. This study significantly provided new insights into the evolutionary patterns and molecular functions of CDP genes, and laid a foundation for further studying the regulatory mechanisms of plant rhizome sprouting.

Origin and Evolution of Marsupial-specific Imprinting Clusters Through Lineage-specific Gene Duplications and Acquisition of Promoter Differential Methylation.

Genomic imprinting is a parent-of-origin-specific expression phenomenon that plays fundamental roles in many biological processes. In animals, imprinting is only observed in therian mammals, with ∼200 imprinted genes known in humans and mice. The imprinting pattern in marsupials has been minimally investigated by examining orthologs to known eutherian imprinted genes. To identify marsupial-specific imprinting in an unbiased way, we performed RNA-seq studies on samples of fetal brain and placenta from the reciprocal cross progeny of two laboratory opossum stocks. We inferred allele-specific expression for >3,000 expressed genes and discovered/validated 13 imprinted genes, including three previously known imprinted genes, Igf2r, Peg10, and H19. We estimate that marsupials imprint ∼60 autosomal genes, which is a much smaller set compared with eutherians. Among the nine novel imprinted genes, three noncoding RNAs have no known homologs in eutherian mammals, while the remaining genes have important functions in pluripotency, transcription regulation, nucleolar homeostasis, and neural differentiation. Methylation analyses at promoter CpG islands revealed differentially methylated regions in five of these marsupial-specific imprinted genes, suggesting that differential methylation is a common mechanism in the epigenetic regulation of marsupial imprinting. Clustering and co-regulation were observed at marsupial imprinting loci Pou5f3-Npdc1 and Nkrfl-Ipncr2, but eutherian-type multi-gene imprinting clusters were not detected. Also differing from eutherian mammals, the brain and placenta imprinting profiles are remarkably similar in opossums, presumably due to the shared origin of these organs from the trophectoderm. Our results contribute to a fuller understanding of the origin, evolution, and mechanisms of genomic imprinting in therian mammals.

Discovery of specific catalytic activity toward IAA/FA by LaSABATHs based on genome-wide phylogenetic and enzymatic analysis of SABATH gene family from Larix kaempferi

The SABATH methyltransferases catalyze methylation of small-molecule metabolites, which participate in plant growth, development and defense response. Given lack of genome-wide studies on gymnosperms SABATH family, the formation and functional differentiation mechanism of the Larix kaempferi SABATH gene family was systematically and exhaustively explored by analyzing gene sequence characteristics, phylogenetic relationship, expression pattern, and enzyme activities. Phylogenetic analysis showed that 247 SABATH genes from 14 land plants were divided into 4 clades, and lineage-specific gene duplication events were important factors that contributed to the evolution of the SABATH gene family in gymnosperms and angiosperms. Substrate specificity analysis of 18 Larix SABATH proteins showed that LaSABATHs could catalyze O-methylation of indole-3-acetic acid (IAA) and farnesic acid (FA), N-methylation of theobromine, and S-methylation of thiobenzoic acid. Furthermore, only LaSABATH2 and LaSABATH29 could catalyze O-methylation of FA, and only LaSABATH30 could catalyze O-methylation of IAA. Homology modeling and molecular docking studies showed the hydrogen bond formed between the His188 of LaSABATH30 and IAA and the noticeable hydrophobic IAA-binding pocket may be helpful for IAA methylation. In this study, identification of proteins with significant specific catalytic activity toward FA and IAA provided high-quality candidate genes for forest genetics and breeding.

Evolution of the nitric oxide synthase family in vertebrates and novel insights in gill development.

Nitric oxide (NO) is an ancestral key signalling molecule essential for life and has enormous versatility in biological systems, including cardiovascular homeostasis, neurotransmission and immunity. Although our knowledge of NO synthases (Nos), the enzymes that synthesize NO in vivo, is substantial, the origin of a large and diversified repertoire of nos gene orthologues in fishes with respect to tetrapods remains a puzzle. The recent identification of nos3 in the ray-finned fish spotted gar, which was considered lost in this lineage, changed this perspective. This finding prompted us to explore nos gene evolution, surveying vertebrate species representing key evolutionary nodes. This study provides noteworthy findings: first, nos2 experienced several lineage-specific gene duplications and losses. Second, nos3 was found to be lost independently in two different teleost lineages, Elopomorpha and Clupeocephala. Third, the expression of at least one nos paralogue in the gills of developing shark, bichir, sturgeon, and gar, but not in lamprey, suggests that nos expression in this organ may have arisen in the last common ancestor of gnathostomes. These results provide a framework for continuing research on nos genes' roles, highlighting subfunctionalization and reciprocal loss of function that occurred in different lineages during vertebrate genome duplications.

Deep evaluation of the evolutionary history of the Heat Shock Factor (HSF) gene family and its expansion pattern in seed plants.

Heat shock factor (HSF) genes are essential in some of the basic developmental pathways in plants. Despite extensive studies on the structure, functional diversification, and evolution of HSF genes, their divergence history and gene duplication pattern remain unknown. To further illustrate the probable divergence patterns in these subfamilies, we analyzed the evolutionary history of HSF genes using phylogenetic reconstruction and genomic syntenic analyses, taking advantage of the increased sampling of genomic data from pteridophytes, gymnosperms and basal angiosperms. We identified a novel clade that includes HSFA2, HSFA6, HSFA7, and HSFA9 with a complex relationship, which is very likely due to orthologous or paralogous genes retained after frequent gene duplication events. We hypothesized that HSFA9 derives from HSFA2 through gene duplication in eudicots at the ancestral state, and then expanded in a lineage-specific way. Our findings indicate that HSFB3 and HSFB5 emerged before the divergence of ancestral angiosperms, but were lost in the most recent common ancestors of monocots. We also presumed that HSFC2 derives from HSFC1 in ancestral monocots. This work proposes that during the radiation of flowering plants, an era during which there was a differentiation of angiosperms, the size of the HSF gene family was also being adjusted with considerable sub- or neo-functionalization. The independent evolution of HSFs in eudicots and monocots, including lineage-specific gene duplication, gave rise to a new gene in ancestral eudicots and monocots, and lineage-specific gene loss in ancestral monocots. Our analyses provide essential insights for studying the evolutionary history of this multigene family.

Evolutionary and Expression Analysis of MOV10 and MOV10L1 Reveals Their Origin, Duplication and Divergence.

MOV10 and MOV10L1 both encode ATP-dependent RNA helicases. In mammals, MOV10 and MOV10L1 participate in various kinds of biological contexts, such as defense of RNA virus invasion, neuron system, germ cell and early development. However, mov10 and mov10l1 in zebrafish are obscure and the evolutionary relationships of mov10 among different species remain unclear. In this study, we found MOV10 and MOV10L1 had some variations despite they possessed the conserved feature of RNA helicase, however, they may originate from a single ancestor although they shared limited homology. A single MOV10L1 gene existed among all species, while MOV10 gene experienced lineage-specific intra-chromosomal gene duplication in several species. Interestingly, the mov10 gene expanded to three in zebrafish, which originating from a duplication by whole genome specific duplication of teleost lineage followed by a specific intra-chromosome tandem duplication. The mov10 and mov10l1 showed distinct expression profiles in early stages, however, in adult zebrafish, three mov10 genes exhibited similar diverse expression patterns in almost all tissues. We also demonstrated mov10 genes were upregulated upon virus challenge, highlighting they had redundant conserved roles in virus infection. These results provide valuable data for the evolution of MOV10 and MOV10L1 and they are important to the further functional exploration.

De Novo Genome Assembly Highlights the Role of Lineage-Specific Gene Duplications in the Evolution of Venom in Fea's Viper (Azemiops feae)

Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea’s Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A2 (PLA2) proteins expressed by genes that likely arose from lineage-specific PLA2 gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families.

Distribution and the evolutionary history of G-protein components in plant and algal lineages.

Heterotrimeric G-protein complexes comprising Gα-, Gβ-, and Gγ-subunits and the regulator of G-protein signaling (RGS) are conserved across most eukaryotic lineages. Signaling pathways mediated by these proteins influence overall growth, development, and physiology. In plants, this protein complex has been characterized primarily from angiosperms with the exception of spreading-leaved earth moss (Physcomitrium patens) and Chara braunii (charophytic algae). Even within angiosperms, specific G-protein components are missing in certain species, whereas unique plant-specific variants-the extra-large Gα (XLGα) and the cysteine-rich Gγ proteins-also exist. The distribution and evolutionary history of G-proteins and their function in nonangiosperm lineages remain mostly unknown. We explored this using the wealth of available sequence data spanning algae to angiosperms representing extant species that diverged approximately 1,500 million years ago, using BLAST, synteny analysis, and custom-built Hidden Markov Model profile searches. We show that a minimal set of components forming the XLGαβγ trimer exists in the entire land plant lineage, but their presence is sporadic in algae. Additionally, individual components have distinct evolutionary histories. The XLGα exhibits many lineage-specific gene duplications, whereas Gα and RGS show several instances of gene loss. Similarly, Gβ remained constant in both number and structure, but Gγ diverged before the emergence of land plants and underwent changes in protein domains, which led to three distinct subtypes. These results highlight the evolutionary oddities and summarize the phyletic patterns of this conserved signaling pathway in plants. They also provide a framework to formulate pertinent questions on plant G-protein signaling within an evolutionary context.

The Visual Opsin Gene Repertoires of Teleost Fishes: Evolution, Ecology, and Function.

Visual opsin genes expressed in the rod and cone photoreceptor cells of the retina are core components of the visual sensory system of vertebrates. Here, we provide an overview of the dynamic evolution of visual opsin genes in the most species-rich group of vertebrates, teleost fishes. The examination of the rich genomic resources now available for this group reveals that fish genomes contain more copies of visual opsin genes than are present in the genomes of amphibians, reptiles, birds, and mammals. The expansion of opsin genes in fishes is due primarily to a combination of ancestral and lineage-specific gene duplications. Following their duplication, the visual opsin genes of fishes repeatedly diversified at the same key spectral-tuning sites, generating arrays of visual pigments sensitive to the ultraviolet to red spectrum of light. Species-specific opsin gene repertoires correlate strongly with underwater light habitats, ecology, and color-based sexual selection.

Comparative analysis of the superoxide dismutase gene family in Cetartiodactyla.

Cetacea, whales, dolphins and porpoises form an order of mammals adapted to aquatic life. Their transition to an aquatic habitat resulted in exceptional protection against cellular insults, including oxidative and osmotic stress. Here, we considered the structure and molecular evolution of the superoxide dismutase (SOD) gene family, which encodes essential enzymes in the mammalian antioxidant system, in the superorder Cetartiodactyla. To this end, we juxtaposed cetaceans and their closest extant relatives (order Artiodactyla). We identified 94 genes in 23 species, of which 70 are bona fide intact genes. Although the SOD gene family is conserved in Cetartiodactyla, lineage-specific gene duplications and deletions were observed. Phylogenetic analyses show thatthe SOD2 subfamily diverged from a clade containing SOD1 and SOD3, suggesting that cytoplasmic, extracellular and mitochondrial SODs have started down independent evolutionary paths. Specific-amino acid changes (e.g. K130N in SOD2) that may enhance ROS elimination were identified in cetaceans. In silico analysis suggests that the core transcription factor repertoire of cetartiodactyl SOD genes may include Sp1, NF-κB, Nrf2 and AHR. Putative transcription factors binding sites responding to hypoxia were (e.g. Suppressor of Hairless; Su(H)) found in the cetaceanSOD1gene. We found significant evidence for positive selection in cetaceans using codon models. Cetaceans with different diving abilities also show divergent evolution of SOD1 and SOD2. Our genome-wide analysis ofSODgenes helps clarify their relationship and evolutionary trajectory and identify putative functional changes in cetaceans.

Recruitment of toxin-like proteins with ancestral venom function supports endoparasitic lifestyles of Myxozoa

Cnidarians are the oldest lineage of venomous animals and use nematocysts to discharge toxins. Whether venom toxins have been recruited to support parasitic lifestyles in the Endocnidozoa (Myxozoa + Polypodium) is, however, unknown. To examine this issue we variously employed transcriptomic, proteomic, associated molecular phylogenies, and localisation studies on representative primitive and derived myxozoans (Malacosporea and Myxosporea, respectively), Polypodium hydriforme, and the free-living staurozoan Calvadosia cruxmelitensis. Our transcriptomics and proteomics analyses provide evidence for expression and translation of venom toxin homologs in myxozoans. Phylogenetic placement of Kunitz type serine protease inhibitors and phospholipase A2 enzymes reveals modification of toxins inherited from ancestral free-living cnidarian toxins, and that venom diversity is reduced in myxozoans concordant with their reduced genome sizes. Various phylogenetic analyses of the Kunitz-type toxin family in Endocnidozoa suggested lineage-specific gene duplications, which offers a possible mechanism for enhancing toxin diversification. Toxin localisation in the malacosporean Buddenbrockia plumatellae substantiates toxin translation and thus illustrates a repurposing of toxin function for endoparasite development and interactions with hosts, rather than for prey capture or defence. Whether myxozoan venom candidates are expressed in transmission stages (e.g. in nematocysts or secretory vesicles) requires further investigation.

The Evolution of Calcification in Reef-Building Corals.

Corals build the structural foundation of coral reefs, one of the most diverse and productive ecosystems on our planet. Although the process of coral calcification that allows corals to build these immense structures has been extensively investigated, we still know little about the evolutionary processes that allowed the soft-bodied ancestor of corals to become the ecosystem builders they are today. Using a combination of phylogenomics, proteomics, and immunohistochemistry, we show that scleractinian corals likely acquired the ability to calcify sometime between ∼308 and ∼265 Ma through a combination of lineage-specific gene duplications and the co-option of existing genes to the calcification process. Our results suggest that coral calcification did not require extensive evolutionary changes, but rather few coral-specific gene duplications and a series of small, gradual optimizations of ancestral proteins and their co-option to the calcification process.

The UDP-Glycosyltransferase Family in Drosophila melanogaster: Nomenclature Update, Gene Expression and Phylogenetic Analysis.

UDP-glycosyltransferases (UGTs) are important conjugation enzymes found in all kingdoms of life, catalyzing a sugar conjugation with small lipophilic compounds and playing a crucial role in detoxification and homeostasis. The UGT gene family is defined by a signature motif in the C-terminal domain where the uridine diphosphate (UDP)-sugar donor binds. UGTs have been identified in a number of insect genomes over the last decade and much progress has been achieved in characterizing their expression patterns and molecular functions. Here, we present an update of the complete repertoire of UGT genes in Drosophila melanogaster and provide a brief overview of the latest research in this model insect. A total of 35 UGT genes are found in the D. melanogaster genome, localized to chromosomes 2 and 3 with a high degree of gene duplications on the chromosome arm 3R. All D. melanogaster UGT genes have now been named in FlyBase according to the unified UGT nomenclature guidelines. A phylogenetic analysis of UGT genes shows lineage-specific gene duplications. Analysis of anatomical and induced gene expression patterns demonstrate that some UGT genes are differentially expressed in various tissues or after environmental treatments. Extended searches of UGT orthologs from 18 additional Drosophila species reveal a diversity of UGT gene numbers and composition. The roles of Drosophila UGTs identified to date are briefly reviewed, and include xenobiotic metabolism, nicotine resistance, olfaction, cold tolerance, sclerotization, pigmentation, and immunity. Together, the updated genomic information and research overview provided herein will aid further research in this developing field.

Genomic signatures of globally enhanced gene duplicate accumulation in the megadiverse higher Diptera fueling intralocus sexual conflict resolution.

Gene duplication is an important source of evolutionary innovation. To explore the relative impact of gene duplication during the diversification of major insect model system lineages, we performed a comparative analysis of lineage-specific gene duplications in the fruit fly Drosophila melanogaster (Diptera: Brachycera), the mosquito Anopheles gambiae (Diptera: Culicomorpha), the red flour beetle Tribolium castaneum (Coleoptera), and the honeybee Apis mellifera (Hymenoptera). Focusing on close to 6,000 insect core gene families containing maximally six paralogs, we detected a conspicuously higher number of lineage-specific duplications in Drosophila (689) compared to Anopheles (315), Tribolium (386), and Apis (223). Based on analyses of sequence divergence, phylogenetic distribution, and gene ontology information, we present evidence that an increased background rate of gene duplicate accumulation played an exceptional role during the diversification of the higher Diptera (Brachycera), in part by providing enriched opportunities for intralocus sexual conflict resolution, which may have boosted speciation rates during the early radiation of the megadiverse brachyceran subclade Schizophora.

Evolution of SHORT VEGETATIVE PHASE (SVP) genes in Rosaceae: Implications of lineage-specific gene duplication events and function diversifications with respect to their roles in processes other than bud dormancy.

MADS-box genes that are homologous to Arabidopsis SHORT VEGETATIVE PHASE (SVP) have been shown to play key roles in the regulation of bud dormancy in perennial species, particularly in the deciduous fruit trees of Rosaceae. However, their evolutionary profiles in Rosaceae have not yet been analyzed systematically. Here, The SVP genes were found to be significantly expanded in Rosaceae when compared with annual species from Brassicaceae. Phylogenetic analysis showed that Rosaceae SVP genes could be classified into five clades, namely, SVP1, SVP2-R1, SVP2-R2, SVP2-R3 and SVP3. The SVP1 clade genes were retained in most of the species, whereas the SVP2-R2 and SVP2-R3 clades were found to be Maleae- and Amygdaleae-specific (Both of the lineages belong to Amygdaloideae), respectively, and SVP2-R1 was Rosoideae-specific in Rosaceae. Furthermore, 10 lineage-specific gene duplication (GD) events (GD1-10) were proposed for the expansion of SVP genes, suggesting that the expansion and divergence of Rosaceae SVP genes were mainly derived by lineage-specific manner during evolution. Moreover, tandem and segmental duplications were the major reasons for the expansion of SVP genes, and interestingly, tandem duplications, a well-known evolutionary feature of SVP genes, were found to be mainly Amygdaloideae-specific. Sequence alignment, selection pressure, and cis-acting element analysis suggested large functional innovations and diversification of SVP genes in different lineages of Rosaceae. Finally, the different growth cycle of Rosa multiflora and their novel expression patterns of RmSVP genes provided new insights into the functional diversification of SVP genes in terms of their roles in processes other than bud dormancy.

Comparative Plastid Genomics of Non-Photosynthetic Chrysophytes: Genome Reduction and Compaction.

Spumella-like heterotrophic chrysophytes are important eukaryotic microorganisms that feed on bacteria in aquatic and soil environments. They are characterized by their lack of pigmentation, naked cell surface, and extremely small size. Although Spumella-like chrysophytes have lost their photosynthetic ability, they still possess a leucoplast and retain a plastid genome. We have sequenced the plastid genomes of three non-photosynthetic chrysophytes, Spumella sp. Baeckdong012018B8, Pedospumella sp. Jangsampo120217C5 and Poteriospumella lacustris Yongseonkyo072317C3, and compared them to the previously sequenced plastid genome of “Spumella” sp. NIES-1846 and photosynthetic chrysophytes. We found the plastid genomes of Spumella-like flagellates to be generally conserved with respect to genome structure and housekeeping gene content. We nevertheless also observed lineage-specific gene rearrangements and duplication of partial gene fragments at the boundary of the inverted repeat and single copy regions. Most gene losses correspond to genes for proteins involved in photosynthesis and carbon fixation, except in the case of petF. The newly sequenced plastid genomes range from ~55.7 kbp to ~62.9 kbp in size and share a core set of 45 protein-coding genes, 3 rRNAs, and 32 to 34 tRNAs. Our results provide insight into the evolutionary history of organelle genomes via genome reduction and gene loss related to loss of photosynthesis in chrysophyte evolution.

The Globin Gene Family in Arthropods: Evolution and Functional Diversity.

Globins are small heme-proteins that reversibly bind oxygen. Their most prominent roles in vertebrates are the transport and storage of O2 for oxidative energy metabolism, but recent research has suggested alternative, non-respiratory globin functions. In the species-rich and ecologically highly diverse taxon of arthropods, the copper-containing hemocyanin is considered the main respiratory protein. However, recent studies have suggested the presence of globin genes and their proteins in arthropod taxa, including model species like Drosophila. To systematically assess the taxonomic distribution, evolution and diversity of globins in arthropods, we systematically searched transcriptome and genome sequence data and found a conserved, widespread occurrence of three globin classes in arthropods: hemoglobin-like (HbL), globin X (GbX), and globin X-like (GbXL) protein lineages. These globin types were previously identified in protostome and deuterostome animals including vertebrates, suggesting their early ancestry in Metazoa. The HbL genes show multiple, lineage-specific gene duplications in all major arthropod clades. Some HbL genes (e.g., Glob2 and 3 of Drosophila) display particularly fast substitution rates, possibly indicating the evolution of novel functions, e.g., in spermatogenesis. In contrast, arthropod GbX and GbXL globin genes show high evolutionary stability: GbXL is represented by a single-copy gene in all arthropod groups except Brachycera, and representatives of the GbX clade are present in all examined taxa except holometabolan insects. GbX and GbXL both show a brain-specific expression. Most arthropod GbX and GbXL proteins, but also some HbL variants, include sequence motifs indicative of potential N-terminal acylation (i.e., N-myristoylation, 3C-palmitoylation). All arthropods except for the brachyceran Diptera harbor at least one such potentially acylated globin copy, confirming the hypothesis of an essential, conserved globin function associated with the cell membrane. In contrast to other animals, the fourth ancient globin lineage, represented by neuroglobin, appears to be absent in arthropods, and the putative arthropod orthologs of the fifth metazoan globin lineage, androglobin, lack a recognizable globin domain. Thus, the remarkable evolutionary stability of some globin variants is contrasted by occasional dynamic gene multiplication or even loss of otherwise strongly conserved globin lineages in arthropod phylogeny.

Evolutionary dynamic analyses on monocot flavonoid 3\u2032-hydroxylase gene family reveal evidence of plant-environment interaction

BackgroundFlavonoid 3′-hydroxlase (F3’H) is an important enzyme in determining the B-ring hydroxylation pattern of flavonoids. In monocots, previous studies indicated the presence of two groups of F3’Hs with different enzyme activities. One F3’H in rice was found to display novel chrysoeriol-specific 5′-hydroxylase activity. However, the evolutionary history of monocot F3’Hs and the molecular basis for the observed catalytic difference remained elusive.ResultsWe performed genome-wide survey of 12 common monocot plants, and identified a total of 44 putative F3’H genes. The results showed that F3’H gene family had underwent volatile lineage-specific gene duplication and gene loss events in monocots. The expansion of F3’H gene family was mainly attributed to dispersed gene duplication. Phylogenetic analyses showed that monocot F3’Hs have evolved into two independent lineages (Class I and Class II) after gene duplication in the common ancestor of monocot plants. Evolutionary dynamics analyses had detected positive natural selection in Class II F3’Hs, acting on 7 specific amino acid sites. Protein modelling showed these selected sites were mainly located in the catalytic cavity of F3’H. Sequence alignment revealed that Class I and Class II F3’Hs displayed amino acid substitutions at two critical sites previously found to be responsible for F3’H and flavonoid 3′5’-hydroxylase (F3’5’H) activities. In addition, transcriptional divergence was also observed for Class I and Class II F3’Hs in four monocot species.ConclusionsWe concluded that monocot F3’Hs have evolved into two independent lineages (Mono_F3’H Class I and Class II), after gene duplication during the common ancestor of monocot plants. The functional divergence of monocot F3’H Class II has been affected by positive natural selection, which acted on specific amino acid sites only. Critical amino acid sites have been identified to have high possibility to affect the substrate specificity of Class II F3’Hs. Our study provided an evolutionary and protein structural explanation to the previously observed chrysoeriol-specific 5′-hydroxylation activity for CYP75B4 in rice, which may also be true for other Class II F3’Hs in monocots. Our study presented clear evidence of plant-environmental interaction at the gene evolutionary level, and would guide future functional characterization of F3’Hs in cereal plants.

Diversified secondary metabolite biosynthesis gene repertoire revealed in symbiotic dinoflagellates

Symbiodiniaceae dinoflagellates possess smaller nuclear genomes than other dinoflagellates and produce structurally specialized, biologically active, secondary metabolites. Till date, little is known about the evolution of secondary metabolism in dinoflagellates as comparative genomic approaches have been hampered by their large genome sizes. Here, we overcome this challenge by combining genomic and metabolomics approaches to investigate how chemical diversity arises in three decoded Symbiodiniaceae genomes (clades A3, B1 and C). Our analyses identify extensive diversification of polyketide synthase and non-ribosomal peptide synthetase genes from two newly decoded genomes of Symbiodinium tridacnidorum (A3) and Cladocopium sp. (C). Phylogenetic analyses indicate that almost all the gene families are derived from lineage-specific gene duplications in all three clades, suggesting divergence for environmental adaptation. Few metabolic pathways are conserved among the three clades and we detect metabolic similarity only in the recently diverged clades, B1 and C. We establish that secondary metabolism protein architecture guides substrate specificity and that gene duplication and domain shuffling have resulted in diversification of secondary metabolism genes.